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We present analytical next-to-leading order results for the 
correlator of baryonic currents at the three loop level with 
one finite mass quark. We obtain the massless and the HQET 
limit of the correlator from the general formula as particular 
cases. We also give explicit expressions for the moments of 
the spectral density. 
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Baryons are an important source of information on 
QCD interactions. Massive baryons containing a heavy 
quark form a rich family of particles with many proper- 
ties experimentally measured The prominent object 
of a QCD analysis of the properties of baryons is the cor- 
relator of two baryonic currents and the spectral density 
associated with it. Calculations beyond the leading or- 
der have not been done for many interesting cases. The 
massless case has been known since long ago [1). The 
mesonic analogue of our baryonic calculation was com- 
pleted some time ago and has subsequently provided a 
rich source of inspiration for many applications in meson 
physics. 

We report on the results of calculating the as correc- 
tions to the correlator of two baryonic currents with one 
finite mass quark and two massless quarks. We give an- 
alytical results and discuss the order of magnitude of the 
as corrections. The massless and HQET limits are ob- 
tained as special cases. We finally present explicit results 
on the moments of the spectral density associated with 
the correlator. The discussion of the impact of our new 
correlator result on baryon phenomenology will be left to 
future publications. 

In the present paper we limit ourselves to the simplest 
case with a three-quark current of the form 



(1) 



which has the quantum numbers of a = 1/2^ baryon. 
^ is a finite mass quark field with mass parameter to, 
u and d are massless quark fields, C is the charge con- 
jugation matrix, e°'"^ is the totally antisymmetric tensor 
and a, 6, c are color indices for the SU (3) color group. 
Other baryonic currents with any given specified quan- 
tum numbers are obtained from the current in Eq. by 
inserting the appropriate Dirac matrices. Such additions 
introduce no principal complication into our method of 
calculation. 




FIG. 1. The calculated (a) two-loop and (b-e) three-loop 
topologies 



The correlator of two baryonic currents is expanded as 
z [ (Tj(x)J(0))e^'?^dx = ^.q-U,iq^) + mU,niq'). (2) 



For reasons of brevity we shall present results only for 
the invariant function Ilmiq^)- The invariant function 
Ilm{q^) can be represented compactly via the dispersion 
relation 



n™(g') = / 



Pmis)ds 



q^ 



(3) 



where Pm{s) is the spectral density. All quantities are un- 
derstood to be appropriately regularized. Since the spec- 
tral density is the real object of interest for phenomeno- 
logical applications, we limit our subsequent discussion 
to the spectral density instead of the correlator. To sim- 
plify formulas we remove a trivial but awkward factor 
containing powers of tt and write 



Pm(s) 
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with a conveniently normalized p{s), 



(4) 



poGs) 1 + ^ln £^ . (5) 
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Here /z is the renormahzation scale parameter, m is a 
pole mass of the heavy quark (see e.g. |^) and as = 
as (yu) . The leading order two-loop contribution shown in 
Fig. 1(a) reads 

po{s) = l + 9z-9z^ - + 6z{l + z)liiz (6) 

with z — vr? js. In the MS-scheme the next-to-leading 
order three-loop contribution is given by 
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pM = 9+— z2_9z3 

^+42z-42z^-^zA In(l-z) 
9 9 y 

„ 154 22 2 58 ,\ , 
2+— z--z^--z3jlnz 

+4 Q + 3z - Sz^ - ^z^^ ln(l - z) In z 

+12z (2 + 3Z+ iz^^ Q In^ z - C(2) 

+4 + 12z + 3z2 - iz^^ Li2(z) 

+24z(l + z) (^Li3(z)-C(3)-iLi2(z)lnz^ (7) 

where Li„(z) are polylogarithms and C("-) is Riemann's 
zeta function. The contributing three-loop diagrams are 
shown in Figs. 1(b) to (e). They have been evaluated 
using advanced algebraic methods for multi-loop calcu- 
lations along the lines decribed in Refs. |^,||. For the 
convenience of the reader we briefly sketch our method 
of integration. All diagrams have first been reduced to 
scalar prototypes. The integrals over massless loops have 
been performed (for recurrent integration where possible) 
and one is left with the basic integral 



{k'^ + m'^)°'{q-k) 



2/3 



(8) 



which is a generalization of the standard object G'(a, (3) of 
the massless calculation. The integral V is known analyt- 
ically and suffices to calculate the diagrams in Fig. 1(c) 
and (d). For the calculation of the diagram shown in 
Fig. 1(e) the basic integral V enters as a subdiagram. 
This subdiagram then is represented in terms of a disper- 
sion integral which makes the whole diagram computable 
in terms of the same V with the argument depending 
on the loop momentum. The final step is a finite range 
(convolution type) integration over this internal momen- 
tum with a spectral density of the basic integral V . The 
reduction to scalar prototypes of the diagram shown in 
Fig. 1(e) leads also to a new irreducible block (i.e. a pro- 
totype not expressible in terms of V) which is related 
to a two-loop master (fish) diagram. The result for this 
diagram is taken from Ref. M. 



The result presented in Eq. (^ is the main result of 
this paper: it represents the full next-to-leading order 
solution. Since the anomalous dimension of the current 
in Eq. (|l) is known to two loop order the result shown 
in Eq. (^) completes the ingredients necessary for an anal- 
ysis of the scalar part of the correlator in Eq. (^ at the 
next-to-leading order level. 

We now turn to the analysis of Eq. (0). Two limiting 
cases of interest are the near-threshold and high energy 
asymptotics. With our result given in Eq. (|^) both limits 
can be taken explicitly. The asymptotic expressions can 
be also obtained in the framework of effective theories 
which can be viewed as special devices for such calcula- 
tions. 

In the high energy (small mass) limit z ^ the cor- 
rection reads 

pi(s) =9 + 83z-47r2z + 21nz 

+50z In z + 12z In^ z - 24zC(3) + O(z^). (9) 

This will lead to the small mass expansion of the spectral 
density in the form 



mp{s) = mjj^{p)p^^--^^-^^{s) 



(10) 
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where pmassicssj-^-j jg -(-j-^g result of calculating the correla- 
tor in the effective theory of massless quarks. The rela- 
tion between the pole (or invariant) mass parameter m 
and the MS mass m-^{^j,) reads 

,2 



In 



(11) 



Note that the massless effective theory cannot reproduce 
the mass singularities (terms like zln(z) in Eq. (|9|)). The 
presence of these singularities is an infrared phenomenon 
and can be parametrized with condensates of local op- 
erators. In our case the first correction in Eq. (|^) 
(and, of course, Eq. (|^)) can be found if the pcrturbative 
value of the heavy quark condensate (5'^') taken from the 
full theory is added The composite operator {^'^) 
should be understood within a mass independent renor- 
mahzation scheme such as the MS-scheme. This value 
(perturbatively, {^'^) ^ Inip^ / m^)) cannot be com- 
puted within the effective theory of massless quarks. It 
represents the proper matching between the pcrturba- 
tive expressions for the correlators of the full (massive) 
and effective (massless) theories. This matching proce- 
dure allows one to restore higher order terms of the mass 
expansion in the full theory from the effective massless 
theory with the mass term treated as a perturbation [ p^ . 
This is justified at high energies. Note that the correc- 
tion to the next order (i.e. of order m^/s) can actually be 
found in this manner because it depends only on one lo- 
cal operator and, therefore, the calculation is technically 
feasible. 
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In the near-threshold hmit E ^ with s = {m + E)"^ 
one exphcitly obtains 



IQE' 
5m 



1 + 
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The invariant function pm suffices to determine the 
complete leading HQET behaviour since one has ^pq + 
iTT-Pm — > m(j>-|-l)p for the leading term. In this region the 
appropriate device to compute the limit of the correlator 
is HQET (see e.g. [|lljl|). Writing 



mp''^'{m,E) = C{m/p,as)^p^^^^iE,fi) 



(13) 



we obtain the known result for p^'^^^{E,p) ||l^ with 
matching coefhcient 0(111/ p, as) |lj]. In this case the 
matching procedure allows one to restore the near- 
threshold limit of the full correlator starting from the 
simpler effective theory near threshold |l^ . 

Note that the higher order corrections in E/m to 
Eq. (^2|) can be easily obtained from the explicit result 
given in Eq. (|^). Indeed, the next-to-leading order cor- 
rection in low energy expansion reads 



Ap'^'{m,E) 
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It is a much more difficult task to obtain this result start- 
ing from HQET. 

We now discuss some quantitative features of the cor- 
rection given in Eq. (^. Of interest is whether the 
two limiting expressions (the massless limit expression 
as given in Eq. (p^ and the HQET limit expression in 
Eqs. ([1^) and (|13D) can be used to characterise the full 
function for all energies. 

In Fig. 1^ we compare components of the baryonic spec- 
tral function in leading and next-to-leading order. Shown 
is the ratio pi{s)/po{s). In the following we shall always 
put /z = TO if it is not explicitly written. One can see that 
a simple interpolation between the two limits can give a 
rather good approximation for the next-to-leading order 
correction in the complete region of s. 
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FIG. 2. The ratio pi(s) / po(s) of the next-to-leading correc- 
tion and the leading order term in dependence of the energy 
square s 



Another informative set of observables are moments of 
the spectral density 



Mn = 



^ = TO«--M„ 



with Af„ dimensionless. We find 



M„.Af(0)^ + ^ In £,1+^, 



where 



and 
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n(n-l)2(n-2)2(n-3) 
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(15) 



(16) 



(17) 



(18) 



The coefficients An are rational numbers. The closed 
form expression for the correction (5„ is long. Instead 
we present explicit expressions for the first several mo- 
ments and the differences between consecutive moments 
in Table |. We see that the difference between consecu- 
tive moments is reasonably small. The absolute value of 
the correction itself in the MS-scheme has no particular 
meaning and can be changed by modifying the subtrac- 
tion scheme. In normalizing to a reference moment at 
n = iV, the others can be expressed through 



A/, 



N 



M 



(0) 
TV 



{l + ^((5„-<5^)}. 



(19) 



One now can find the actual magnitude of the correction 
from Tabic ^. Indeed, for any given precision and range of 
n the set of perturbatively commensurate moments can 
be found. 

Note that moments represent massive vacuum bubbles, 
i.e. diagrams without external momenta with massive 
lines. These diagrams have been comprehensively ana- 
lyzed in Refs. |l6|,p7|. The analytical results for the first 
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few moments at three-loop level can be checked indepen- 
dently with existing computer programs (see e.g. ||T^). 

Moments of the spectral density are convenient ob- 
servables for phenomenological applications. Exact re- 
sults for correlators with massive particles are not known 
in higher orders for many important physical channels. 
Therefore one generally uses approximate procedures. 
We consider the two formal limiting cases as a base for 
possible interpolation. Both are simpler than the full 
calculation. 

For the first approximation, the high energy or mass- 
less approximation, we first check the numerical differ- 
ence between the exact result and the massless approx- 
imation. The massless limit for the moments means 
that one formally integrates the expression pmassiossj-g-j 
in Eq. ( p!o[) over the range {m^, oo) to obtain 



A ^massless 

The result is 



oo ^massless 



^6 — 2n^^masslcss 



M, 



massless 

1 

1 

71 — 3 



(20) 
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This is the extrapolation from the side of high energies. 
The opposite limit for the moments is obtained by us- 
ing the near-threshold spectral density for the integration 
along the whole s-axis. This approximation is an extrap- 
olation from the threshold region which is clearly a poor 
approximation far from threshold. However, it formally 
exists for sufficiently large n and can be obtained with 
the simple HQET result shown in Eqs. ( p^ ) and (13). 
The moments are given by 



= 2m^ 



oo th 



oo ^thr 



(m + £;)2"-i ■ 



(22) 



From Eq. (|12| ) we have 

26(2n-8)! 



^thr(O) ^ 
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(2n-2)! 



(l + 0(n-i)) (23) 



where the large n result specifies the region of applicabil- 
ity of the near-threshold approximation for the moments 
of the spectral density. The correction is 



Mf^ =Mf'^(°)|l-|-^(^ln(^. 



(24) 
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41n2-4(7/;(2n-7) - ipi6)) 



where 'il'{z) is Euler's 7/>-function. 

One can see that the corrections to the moments basi- 
cally reflect the shape of the correction to the spectrum 



as given in Eq. (|^). The massless approximation is rea- 
sonably good for relative corrections for the first few mo- 
ments despite the unfavorable shape of the weight func- 
tion 1/s". It can be improved by changing the subtrac- 
tion point p or by resumming the integrand [|l^ which 
lies beyond the scope of finite order perturbation theory 
though. 

To conclude, we have computed the next-to-leading 
perturbative corrections to the finite mass baryon cor- 
relator at three- loop order. Technically, the method al- 
lows one to obtain analytical results for two-point correla- 
tors of composite operators with one finite mass particle 
which can be compared to HQET results. Corrections 
in E/m near threshold are easily available from our ex- 
plicit results. From threshold to high energies the exact 
spectral density interpolates nicely between the leading 
order HQET result close to threshold and the asymptotic 
mass zero result. Going even one order higher it is very 
likely that the full four-loop spectral density can be well 
approximated by the corresponding massless four-loop 
result which can be calculated using existing computa- 
tional algorithms |2^J2ltl. 
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n 




Sn — S„- 


4 
5 


3 

13/2 


3.500000 


6 


17/2 


2.000000 


7 


535/54 


1.407407 


8 


1187/108 


1.083333 


9 


64093/5400 


0.878333 


10 


22691/1800 


0.737037 


11 


1167767/88200 


0.633878 


12 


2433499/176400 


0.555357 


13 


68055703/4762800 


0.493665 


14 


14034047/952560 


0.443968 


15 


348916495/23051952 


0.403114 


16 


8935543717/576298800 


0.368960 


17 


3086442535481/194788994400 


0.340003 


18 


3147830736641/194788994400 


0.315152 


19 


3205021446217/194788994400 


0.293603 


20 


6517078055669/389577988800 


0.274746 


21 


382499185005589/22517607752640 


0.258112 



TABLE I. Values for the rational part An of the first mo- 
ments Sn and their relative difference Sn — Sn-i 
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